Absence of intermediate spin state in hole-doped LaCoOa: dynamical mean-field study 
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The hole and electron doped LaCoOg is studied using the dynamical mean-field theory. The one- 
particle spectra are analyzed and compared to the available experimental data, in particular the 
x-ray absorption spectra. Analyzing the temporal spin-spin correlation functions we find the atomic 
intermediate spin state is not important for the observed Curie- Weiss susceptibility. Contrary to 
commonly held view about the roles played by the t2g and eg electrons we find narrow quasi-particle 
bands of t2g character crossing the Fermi level accompanied by strongly damped Cg excitations. 

PACS numbers: 75.30.Wx,71.28.-|-d,75.20.Hr,71.10.Fd 



Correlations arising from the electron-electron inter- 
action within the narrow d bands put transition metal 
oxides among the most puzzling materials, properties of 
which may change substantially upon a small variation 
of temperature, pressure, carrier concentration or some 
other control parameter. LaCoOa is a prominent exam- 
ple. A small gap non-magnetic insulator at low temper- 
atures, its band structure is well described by an effec- 
tive non-interacting electron pictureii Nevertheless, in- 
creased temperature or hole doping reveal its strongly 
correlated nature. Between 50 and 100 K the magnetic 
susceptibility sharply increases and acquires Curie- Weiss 
(CW) character, signaling the appearance of local mag- 
netic moments.— The charge gap continuously disappears 
between 450 and 600 K accompanied by transfer of the 
optical spectral weight on several eV scaled The gen- 
erally accepted explanation is essentially atomic physics 
of the low spin (i^gCg) ground state with a nearby ex- 
cited state of either high spin (HS, ^Ig^g) or intermedi- 
ate spin (IS, i^gSg) character, where the precise nature 
of the excited state is not settled. The hole doping of 
the Co d bands due to substitution of La with Sr, Ca or 
Ba also leads to a strong magnetic rcsponse^i^ At < 18% 
Sr concentration microscopically inhomogeneous phase is 
observed that can be described as magnetic clusters sep- 
arated by non-magnetic matrix characterized by insulat- 
ing spin-glass behavior at low temperatures. Above 18% 
Sr concentrations the material becomes a homogeneous 
ferromagnetic (FM) metal. The transition between these 
two regimes, generally described as a percolation of the 
magnetic clusters, is associated with insulator-to- metal 
transition. 

Properties of Lai_a;Sra;Co03 are commonly discussed 
in the context of double-exchange model which provides a 
satisfactory description of related Lai-^^Sr^^MnOs family. 
In this picture the t2g electrons are localized on the metal 
atom forming the local spin moment while eg electrons 
form dispersive bands. Numerous comparative studies, 
however, found sizeable differences between cobaltitcs 
and manganites. Colossal magnetoresistance, the hall- 



mark of manganite physics, is not found in cobaltites4i^ 
The NMR relaxation rates in cobaltites are several orders 
of magnitude larger than in manganites.^ The linear spe- 
cific heat coefficient in Lao.TSro.aCoOs is 16 times larger 
than in Lao.ySro.sMnOa.— These observations raise the 
question of the relevance of the double-exchange picture 
for doped cobaltites. 

The strong T-dependence of physical properties even 
in the parent compound LaCoOs make theoretical de- 
scription of cobaltites challenging. Several density func- 
tional and Hartree-Fock calculations on doped cobaltites 
were reported,—^— which cannot treat dynamical elec- 
tronic correlations. In this Letter we use the combina- 
tion of the dynamical mean-field theoryi^ and the den- 
sity functional theory (DFT+DMFT)i^^ to study the 
one-particle spectra and magnetic properties of doped 
LaCoOa. Our results, which are relevant for doping 
X > 0.2 with metallic FM ground state, contradict 
the double-exchange picture and the important role at- 
tributed to the IS atomic states. Our DMFT solution 
contains strongly damped local Cg excitations and dis- 
persive t2g bands crossing the Fermi level. We find that 
the local moments in doped cobaltites originate mostly in 
the HS state and give this statement quantitative mean- 
ing. We will argue that p — c? covalency is necessary to 
explain the observed behavior. 

The calculation proceeds in several steps. First, the 
LDA band structure is determined using WIEN2ki^ and 
an effective Hamiltonian spanning the Co d and O p 
bands is constructedJ^iii For all dopings we used the 
same Hamiltonian obtained for the 5 K lattice parame- 
ters of the stoichiometric LaCoOa. Adding the explicit 
electron-electron interaction within the Co d shell we ar- 
rive at the multi-band Hubbard Hamiltonian of the form 



Pk 



+ }^iyf''. (1) 



Here, (pj^.) is an operator-valued vector whose ele- 
ments are Fourier transforms of dia (Pij), which anni- 
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hilate the Co d (O p) electron in the orbital a (7) in 
the ith unit cell. The on-site interaction Wf''' is ap- 
proximated by the density-density form for the sake of 
computational efficiency. The results of our previous 
study on SrCoO^i^ justify this approximation. The in- 
teraction parameters U=6 eV and ,/ =0.8 eV were taken 
from Ref. ITol . The double-counting term E^c approx- 
imately corrects for the explicitly unknown mean-field 
part of the interaction coming from LDA. It was chosen 
to equal the orbitally averaged Hartree part of the self- 
energyi^ The hybridization expansion continuous-time 
quantum Monte Carlo^ii^ with the improved estimator 
for the self-energ}*2^ was used to solve the auxiliary im- 
purity problem. Since the perturbation order is strongly 
imbalanced between weakly hybridized t2g electrons and 
strongly hybridized eg electrons, we use an improved 
measurement based on continuous-estimators to avoid 
poor statistics of the t2g Green's function. All the real- 
axis spectra were obtained by direct analytic continua- 
tion of the self-energy using maximum entropy methodi^ 
The calculations were performed at the temperature of 
580 K. 

We start the presentation of our results with the anal- 
ysis of local quantities in the non-spin-polarized state. 
The undoped system has the one-particle spectrum of 
a band insulator with a crystal-field gap between the 
filled t2g and empty Cg statesii^ The p d covalency 
leads to actual occupancies of the t2g and Cg orbitals of 
5.91 and 0.84, respectively. Fig. [T] shows the variation 
of the orbital occupancies with doping. The strong non- 
monotonicity indicates a departure from an independent 
particles behavior, in particular on the hole-doped side, 
where approximately two Cg electrons and three t2g holes 
are created for each hole added to the system. This is 
accompanied by growing CW spin response as reflected 
in the local spin-spin correlation functions (right panel 
of Fig. [T]). The competition between the different spin 
states of Co atom in LaCoOa is sensitive to the rela- 
tive positions of the Co d and O p bands as the p — Cg 
hybridization contributes substantially to the apparent 
crystal-field splitting. It is therefore crucial to establish 
that the present results do not depend on the variation 
of the double-counting correction E^c, the definition of 
which is somewhat fuzzy in any DFT-|- approach. In ad- 
dition to the choice of self-consistently adjusted E^c, used 
throughout the paper (sec inset in Fig. [Ij, we have per- 
formed a scries of test calculations with E^c fixed to the 
value taken from the undoped system (stars in Fig. [T]). 
Although the two values of E^c differ as much as 1.3 eV 
at x=0.3 we observe still the same behavior of the orbital 
occupancies. 

Associating the CW response of an atom embedded in 
solid with a particular atomic multiplet may not be pos- 
sible in a strongly hybridized and even metallic system. 
In Ref. [1^ we have introduced the atomic state correla- 
tion matrix 11^^ = dr(Pa(T)P^(0)) to address this 
question in undoped LaCoOs. IIq,^ includes effects of 
both quantum and statistical fluctuations and provides 
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FIG. 1: (color online) Left panel: the doping dependence 
(hole doping < 0) of the orbital occupancy. The stars mark 
results obtained with fixed double-counting correction. Right 
panel: the local spin-spin correlation functions: black (full 
line) - the undoped system, blue (dash-dotted line) - hole 
doping x=0.1, 0.2, 0.3 (corresponding to increasing magni- 
tude), and green (dashed line) - electron doping of 0.2. The 
inset shows the doping dependence of _Edc. 



LS 

LS+eg 
LS+2eg 

IS 

HS 

HS+e„ 



LS+en LS+2e„ 



HS HS+eg HS+tjg 



HS+t 



■2g 



m 


)06 


0.004 


0.014 






0.044 


0.111 






1.751 



FIG. 2: Atomic state correlation matrix Yla/i for the hole dop- 
ing of x=0.2 with the contribution of the most abundant d^, 
d'' and states (see Ref. [l^ for the notation). The num- 
bers show the contributions to the local susceptibility Txioc 
summed over the blocks. 



information about the likelihood of finding the Co atom 
in state a (column sum) as well as about the average 
duration of the visit to a given state (large diagonal ele- 
ments mean long visits). Weighted with the spin matrix 
elements, S'^S'^IIq,^ is the contribution of the pair a/? to 
the local susceptibility. The correlation matrix in Fig. [2] 
reveals that hole doping causes increasing weight of the 
HS block which dominates the local susceptibility, sim- 
ilar to the thermal effect in stoiciometric LaCoOsj ^^i^^ 
This observation is consistent with our previous study 
on SrCoOa-i^ 

The one-particle spectral functions decomposed into 
their momentum and orbital contributions are shown in 
Fig. [21 The stoichiometric system exhibits sharp bands 
with well defined dispersion, rcfiecting the fact that an 
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FIG. 3: (color online) Left panel: orbitally resolved spectral 
density: Co eg (black), Co t2g (red), O p (shaded blue), p — n 
(lighter blue) for various dopings. Right panel: the corre- 
sponding k-resolved spectral functions along the high sym- 
metry lines in the Brillouin zone depicted as color plots. 



additional eg electron or hole behaves as a free par- 
ticle. The electron doping introduces also some corre- 
lation effects, but we focus on the hole doping which is 
relevant for Lai_a;Sra;Co03. At x=0.1 the chemical po- 
tential is pinned close to the top of the t2g band (not 
shown). The Cg spectral weight is transferred from the 
coherent Cg — p anti-bonding band into incoherent struc- 
tures around the chemical potential and fills the gap. 
With the number of holes growing the chemical poten- 
tial shifts deeper into the t2g band, which exhibits in- 
creasing mass renormalization. At higher doping lev- 
els some Cg spectral weight builds up at the chemical 
potential. These are strongly damped excitations with 
a lifetime 15 times shorter than that of t2g excitations 
(-ImEe,(0) ~ 0.6 eV, -ImEt,^(0) « 0.04 eV). For the 
full evolution of the hole doped spectra up to a:: = 0.3 see 
the Supplemental Materiali^ 

The overall shape of the O p spectrum is unaffected by 
doping. In particular, the positions of the non-bonding 
(^ —3 eV) and the cr-bonding peaks —5.5 eV) remain 
unchanged. This compares well with the photoemission 
data of Saitoh et al.^ where the peaks shift as little as 
0.4 eV between x = and x = 0.4, and a posteriori jus- 
tifies our choice of the double-counting correction. The 
calculated as well as the experimental electron-removal 
(photoemission) spectra exhibit only a minor variation 
with the hole doping. The electron-addition spectrum 
experiences more pronounced changes. These can be 
measured by means of the x-ray absorption spectroscopy. 
Small core-hole effects in the oxygen K-edge spectrum al- 
low a direct comparison of the XAS data with the O p 
spectral function. In Fig. 2] we show the evolution of 
the positive energy part of the O p spectra broadened 
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FIG. 4: Left: the unoccupied O p spectra (broadened with 
Lorentzian) for various hole dopings. Right: The experimen- 
tal XAS spectra reproduced from Ref. [l^. 



to account for XAS linewidth. The observed transfer 
of spectral weight is driven by the corresponding evo- 
lution of the Co eg spectrum. Our spectra provide a 
good description of the recent experimental observations 
of Mcdling et alw^ In the LSDA theory used by Medling 
et al. the hole doping shifts the chemical potential deeper 
into the t2g peak while the bandstructure remains essen- 
tially rigid. As a result some O p states are brought above 
the chemical potential, due to hybridization with the Co 
t2g's, leading to the increase of the lower energy peak. 
The peak at higher energy, coming from eg—p hybridiza- 
tion is unaffected by doping in the theory of Medling et 
al. This is at odds with the experimental data, which 
show an increase of the low-energy peak with a simulta- 
neous decrease of the high-energy one, accompanied by 
appearance of an isosbestic point<^ The spectral densi- 
ties obtained with the present theory in Fig. |3] exhibit 
precisely this behavior. The increase of the low-energy 
peak in our data arises partly from the depopulation of 
the t2g states and partly from the incoherent eg excita- 
tions, which gradually fill the gap in LaCoOa spectrum. 
Decrease of the high-energy peak in Fig. 2] follows the 
diminishing weight of of the eg anti-bonding peak. 

For the hole doping x = 0.3 there is a stable FM so- 
lution at T=580 K. The ordered moment of 1.2 fiB per 
formula unit is close to saturation as was confirmed by 
calculation at lower temperature. The spectral functions 
are shown in Fig. [Sj There is a clear similarity between 
this spectrum and the spectrum calculated for SrCoOa 
In particular, the majority t2g and minority eg states are 
absent from the chemical potential, while the minority t2g 
states form a well defined band, which is the main con- 
duction channel. The majority eg states form a strongly 
damped band which crosses the chemical potential. The 
spin-averaged spectra exhibit only minor difference to the 
PM spectra, while the eg lifetime is by a factor of two 
larger in the FM phasei^ The atomic state correlation 
matrices show no sizable redistribution of the multiplet 
weights, in particular no enhancement of the IS block, 
due to the spin polarization. This shows that inter-site 
FM correlations, which are completely absent in the PM 
solution, do not stabilize the IS state. 

Lai_a;Sra;Co03 is commonly compared to perovskite 
manganites and the double-exchange model is invoked 



4 




FIG. 5: (color online) Left panel: orbitally resolved spectral 
density: Co Eg (black), Co t2g (red), O p (shaded blue), p — 
TT (lighter blue) for the ferromagnetic solution obtained at 
r=580 K and hole doping x=0.3. (Minority spin densities 
are multiplied with —1). Right panel: The corresponding It- 
resolved spectra for the majority (top) and minority (bottom) 
polarization. 



with the partially populated eg band providing the con- 
duction channel. In this setting the if IS state is a nat- 
ural building block maximizing the kinetic energy gain. 
Our results provide a picture, which differs in two im- 
portant aspects: (i) the HS state is identified as the ma- 
jor contributor to the CW susceptibility, (ii) t2g quasi- 
particles are found at the Fermi level while the Cg exci- 
tations are strongly damped. In the rest of the paper we 
discuss our results in the light of the known experimental 
data. 

First, the computed saturated FM moment at a; = 0.3 
is 1.2 hb per formula unit, slightly below the experimen- 
tal valuej^i^ although our analysis shows that it origi- 
nates largely from the HS state, i.e. the saturation mag- 
netization is not a good indicator for the underlying mul- 
tiplet dynamics. Second, we studied only homogeneous 
solutions with full crystal periodicity. The results are, 
therefore, not relevant for the spin-state polaron regime 
at very low dopingi^SiSi It is plausible that while IS states 
participate in the polaron formation they are not impor- 
tant at higher doping levels. The polaron is stable if the 
kinetic energy gain from delocalizing an Cg electron over 
the polaron volume outweighs the cost of exciting the 
participating atoms into the IS state. While the latter 
term grows linearly with the size of the polaron, the ki- 
netic energy gain is largest for small polarons and quickly 
saturates as the polaron grows. Therefore only small po- 
larons are stable and a substantial change of magnetic 
properties around x = 0.05^ suggests that quite differ- 
ent microscopic physics is involved. 

The other arguments against the IS state in the metal- 
lic phase (x > 0.2) are the lack of substantial magne- 
toresistance (found in manganites) ,— no Griffiths phase,— 
and observation of positive curvature of inverse suscepti- 
bility above The last observation points to increase 
of the paramagnetic moment with temperature, while in 



the IS scenario one would expect decrease of IS atoms 
in magnetically disordered phase. Instead of the double- 
exchange mechanism stabilizing the spin-state polarons, 
which relies on the spin polarization on nano-scale, we 
attribute the population of HS state to the following co- 
valency effect. Introducing a hole on Co atom is likened 
to substitution of a more electronegative atom in the Co 
lattice, which is known to induce the HS state on its 
neighbors2ii2^ due to reduction of the hybridization part 
of the crystal field. Self-consistent adjustment of the sys- 
tem to this situation establishes a new balance of t2g 
holes and Cg electrons. 

The second, and perhaps most important, result of the 
present study is the observation of much stronger damp- 
ing of Eg excitations than their t2g counterparts. To pro- 
vide a qualitative explanation of this effect we make the 
following observations. First, the effect is still present in 
spin-polarized solution, which is easier to analyze since 
the occupied majority t2g and empty minority Cg orbitals 
do not play active role. Second, we found a similar ef- 
fect in a simplified two-band model of Ref. [s^ when hole 
doped. The situation can be viewed as two partially filled 
asymmetric bands of spinless fermions. The second Born 
approximation for the self-energy in the imaginary time 

I],(t) = (72G<,(t)G_,(t)C(/3-t), (2) 

where a = ±1 indexes the two bands, shows that the 
ratio of the self-energies at /3/2 is the inverse ratio of 
the Green functions. This approximately means that 
the band with lower density at the Fermi level has more 
strongly damped quasiparticles. In other words, the slow 
t2g holes act as 'static' scattering centers for the fast Cg 
electrons, while the t2g holes perceive the Cg electrons as 
'homogeneous gas'. Although calculation of the trans- 
port coefficients is beyond the scope of this study we 
make several speculative comments. The present solu- 
tion is consistent with the absence of colossal magne- 
toresistance. The anomalous Hall effect observed in FM 
metallic phaseS^ may be related to the spin-orbit effects 
in the t2g quasiparticle band crossing the Fermi level. 
Finally, the contribution of both t2g and eg channels 
with different quasiparticle dampings, of which the Cg 
one is more sensitive to the magnetic order, should be 
considered when interpreting the T-dependence of ther- 
mopowerisi 

In summary, we have used dynamical mean-field theory 
to investigate hole doped LaCoOs. Our results provide 
a good description of the spectral and magnetic prop- 
erties of the homogeneous metallic phase for x > 0.2. 
We have found that IS state does not play important 
role in this regime and that the CW susceptibility can 
be attributed to HS state and fiuctuation around it. We 
found coherent t2g band crossing the Fermi level, while 
the Cg excitations are strongly damped and appear at the 
Fermi level only at dopings x > 0.2. This leads to the 
conclusion that double-exchange model with Cg band as 
the only conduction channel is not appropriate for doped 
cobaltitcs. 
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